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1. Introduction 

The gauge gravity correspondence, |T], provides a fascinating tool to explore strongly 
coupled field theory Recently, it has been applied to condensed matter systems 
yielding interesting qualitative results for systems exhibiting a superconducting phase 
[[J. In these, the bulk "classical" theory has gauge and charged scalar (or fermi) 
fields, and a black hole provides a finite temperature. Typically, no hair theorems, 
[[|, would lead us to believe that the scalar field must be in its vacuum even in a 
highly charged black hole background, however, such theorems do not take account of 
the fact that in anti de Sitter spacetime scalar fields can have an apparently tachyonic 
mass, provided it is not too large, H|; the confining properties of anti-de Sitter (adS) 
spacetime in essence prevent the instability from setting in unless the wavelength 
is sufficiently small. Once the scalar has a negative mode, the usual conditions of 
no-hair theorems cease to hold, and it becomes possible for the scalar to condense 
out of the vacuum in an analogous fashion to the condensation of SU(2) t'Hooft 
Polyakov fields outside a magnetic Reissner- Nordstrom solution |5|] . 

The basic picture is that at a sufficiently low temperature, it becomes ener- 
getically feasible for a charged scalar to acquire an expectation value near the event 
horizon of the black hole, 0, spontaneously breaking the gauge symmetry and screen- 
ing the charge and mass of the black hole. Because the asymptotic true vacuum is 
symmetric, the scalar has a power law fall-off near the boundary, and the coefficient 
of this fall-off can be interpreted as a condensate in the boundary theory. There is a 
reasonably complete understanding of these systems, including varying scalar mass 
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and potential, the gauge group, the number of spacetime dimensions, as well as hav- 
ing magnetic fields present and the stability of the system 0, [7], |], [| [HJ [11], [12], [L3|] . 
While most of these models are "bottom up" , in the sense of being empirically con- 
structed, similar models have been found from the top down perspective [nj, and it 



is therefore of interest to consider more general stringy aspects of these models. 
In a previous paper, some of us, JT5| , explored the stability of these empirical 



models to leading order corrections, by including on the gravitational side higher 



curvature terms, specifically, the Gauss-Bonnet (GB) invariant |16| which is believed 
to be the 0(a') correction to low energy string gravity (17]]. We found that many of 
the qualitative features of the holographic superconductors were stable under higher 
curvature corrections, however, upon studying the conductivity, a more interesting 
story emerged. In || , Horowitz and Roberts studied 2 and 3 + 1-dimensional holo- 
graphic superconductors for a variety of bulk scalar masses. They discovered that 
the energy gap typically present in the real part of the conductivity, u g , was always 
about eight times as large as the critical temperature of the superconductor: 

co 9 ~ 8T C (1.1) 

this suggested that in spite of the different dimensions of the boundary condensate, 
there was some universal mechanism governing their superconductivity. Upon adding 
in higher curvature corrections however, this "universal" relation completely disap- 
peared. This analysis was confirmed for more general models in |18[, and although 



the expressions for conductivity used in these papers are not strictly accurate (see 
comments in section |4] and appendix [AJ) the gap result is. The analyses presented in 
JT5| , [T8| are however in the probe limit, i.e. where the matter fields do not backreact 
gravitationally on the spacetime, and since backreaction does typically alter the crit- 
ical temperature, it is likely that a full computation will alter the relation between 
Ug and T c , potentially restoring the universal relation, at least for some value of 
the gravitational coupling for each a. In [19], the stability of the scalar condensate 



to higher curvature corrections was tested, and it was found that backreaction low- 
ered T c , as with Einstein gravity, thus making condensation harder. There was no 
qualitative difference however with the probe results. 

In this paper we include the full effect of gravitational backreaction on the holo- 
graphic superconductor with higher curvature terms. We find analytic bounds on the 
critical temperature, and cross-check these against exact results obtained by numer- 
ical computation. We then numerically compute the conductivity, and demonstrate 
that the conductivity does not in fact have a universal gap. Indeed, we find that even 
in Einstein gravity there is no universal gap once backreaction is taken into account. 
The organization of the paper is as follows: We first review the GB holographic su- 
perconductor in section 0, discussing the bulk model, clarifying possible ambiguities 
in the choice of mass of the scalar field and setting out some of the basic properties of 
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solutions to the system of equations. In section ^|we analyse the bulk equations ana- 
lytically, providing bounds on the critical temperature and present numerical results 
for the condensation and critical temperature of the holographic superconductor. We 
then investigate the conductivity in section [|, first deriving the backreacted conduc- 
tivity equation for GB gravity, then solving it numerically. Finally, we conclude in 
section 



2. The Gauss-Bonnet bulk superconductor 

We begin with the Einstein-Gauss-Bonnet (EGB) gravitational action coupled to a 
massive charged complex scalar field and a U(l) gauge field: 
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where g is the determinant of the metric, and R a bcd, Rab and R are the Riemann 
curvature tensor, Ricci tensor, and the Ricci scalar, respectively. We take the Gauss- 
Bonnet coupling constant a to be positive, and the negative cosmological constant 
term, —6/L 2 , has been written in terms of a length scale, L. Note k 2 = 8ttG^, gives 
an explicit Planck scale, and q is the charge, and m the mass, of the scalar field ip 1 . 
The equations of motion can be readily derived as: 
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T ab is the matter energy momentum tensor 
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(2.3) 

(2.4) 



and T> a = V a — iqA a is the gauge covariant derivative. 

In the absence of any matter sources, the EGB equations have a pure adS solution 



ds z 
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1 Note that we follow Horowitz et al., [Q, in taking a purely quadratic potential for the scalar 
field. 
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Thus the actual curvature of the adS spacetime is renormalized away from the cos- 
mological constant scale, L, once a is nonzero. Since L e < L, one could interpret 
switching a on as strengthening gravity, in the sense that a shorter lengthscale cor- 
responds to stronger curvature. 

To examine holographic superconductivity, we look for plane-symmetric black 
hole solutions with or without a nontrivial scalar, but with a nonzero charge. Taking 
the metric ansatz 

ds 2 = f(r)e 2 ^dt 2 - ^ - ^Jdx 2 + dy 2 + dz 2 ) (2.7) 

f{r) L 2 e 

in the absence of a scalar field there is an analytic charged black hole solution, [^(J . 
with v = 0, and 

A = cf>(r)dt = |- (l - ^ dt (2.8) 

where Q is the charge of the black hole (up to a geometrical factor of 4ir), and 
r + is the event horizon, which determines the "ADM" mass of the black hole [pT[|. 
In order to avoid a naked singularity, we need to restrict the parameter range as 
a < L 2 /4. Note that in the Einstein limit (a — > 0), the solution ( |2.9| ) goes to the 
Reissner-Nordstrom adS black hole: 



rl(r^_r_l\ 2k 2 Q 2 (r% r 2 



f(T) = r± -T-^ +^lM4--± • (2-10) 



In the Chern-Simons limit, a = L 2 /4, (so called because in odd dimensions the 
gravitational lagrangian becomes the potential for the Euler density in one dimension 



up, see e.g. P2| ) the Newtonian potential takes the simpler form: 




/M = ^- T fJl + 2L^(l-^). (2.11) 



The superconducting phase corresponds to a "hairy" black hole, where the scalar 
field has condensed out of its symmetric state and screens the charge of the black 
hole. We can see that this will happen at sufficiently low temperature from Gubser's 
rough argument, ||, using the scalar 'effective mass' m 2 s = m 2 — q 2 <ft(r) 2 / f(r). For 
low temperature black holes, f(r) increases slowly away from the horizon, giving a 
large negative mass squared over a sufficient range for an instability to set in. This 
is of course confirmed by the numerical results. 

For the purposes of our investigation, we wish to have a fixed mass in order 
to focus on the effects of backreaction and the higher curvature terms. In |7J, the 



- 4 - 



mass of the scalar was given in terms of the adS lengthscale (m 2 = —(D — 2)/L 2 
in .D-dimensions), and thus fixed with respect to the cosmological constant via the 
Einstein equations. In |L5], |19[ this same value of the mass was also chosen with 
respect to the cosmological constant scale L. However, in EGB gravity, the adS 
lengthscale is dependent on both L and a via ( |2.6| ), and fixing the mass with respect 
to L and not L e means that the dimension of the operator in the dual boundary 
theory varies with a. It is possible therefore that some of the phenomena observed 
in [[15|, |19|] are a consequence of this varying dimension, rather than intrinsic to the 
system, and we therefore fix the mass relative to the asymptotic adS lengthscale, 
m 2 = — 3/Lg, in order that the boundary operator has fixed dimension 3. 

In order to look for the hairy black hole solution, we take the standard static 
ansatz for the fields: 

A a = <f>(r)8° a , = (2.12) 

where without loss of generality if) can be taken to be real. The full system of gravity 
and gauge-scalar equations of motion is then obtained as 2 : 

0"+ 0-^-2^0 = 0, (2.13) 

1 — r] v = ~^~ r ^ + ( 2 - 15 ) 
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where a prime denotes derivative with respect to r. These equations have several 
scaling symmetries, similar to those noted in J7J, although as we have explicitly kept 
the Planck scale, an additional symmetry corresponding to a rescaling of energy is 
present. 

1. r — > ar, t,x i — > at, ax 1 , L — > aL, q — >■ q/ 'a, a — >■ a 2 a, A — > aA. 

2. r -^br,t -> t/b, x i -> x i /b, f -> b 2 f, ->■ b<p. 

3. — > c(f>, if) — > cif), q — y q/c, k 2 — y k 2 /c 2 . 

We use these rescalings to set L = Q = q = 1 for numerical convenience. Note that 
in [0, the Planck scale was set to unity a priori, hence the probe limit corresponded 
to q — y oo, as can be seen from the third rescaling. Also in contrast to pi, we choose 



2 Note, the ij component of the EGB equations is not independent via a Bianchi identity. 
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to fix Q = 1, so that we are explicitly holding the charge parameter fixed in all our 
computations. 

Finally, the Hawking temperature is given by 



T = i/'(r)e"« 



(2.17) 



r=rj r 



this will be interpreted as the temperature of the conformal field theory (CFT). 

In order to solve our equations we need to impose boundary conditions at the 
horizon and the adS boundary. 

• Horizon: 

The position of the horizon, r + , is defined by f(r + ) = 0. Demanding regularity 
of the matter fields and metric at the horizon gives: 

2 

Tfl 

0(r + ) = 0, V'(r+) = 777^-)^( r +) " ( 2 ' 18 ) 

Then equations (|2.15|) and ( |2.16| ) give: 

,/ s 2/t 2 / „ , 2 0'(r+) 2 ^(r + ) 2 \ 

v (r+) = - r+ V (r+) + f \r!yXti) ) ^ 

• Boundary: 

As we want the spacetime to asymptote to adS in standard coordinates, we look 
for a solution with 

v -)■ , f(r) ~ — t as r — >• oo . ( 2 -21) 
Asymptotically the solutions of and ip are then found to be: 

<Kr)~P~%> ^(r)~^- + ^, (2.22) 

where A-t = 2 ± v/4 + m 2 L^ for a general mass m 2 . In order to have a normal- 
izable solution we take C_ = 0; P and C + are then fixed by consistency with the 
near horizon solution. According to the adS/CFT correspondence, we can interpret 
(C A+ ) = C + , where Oa + is the operator with the conformal dimension A + dual 
to the scalar field. As already mentioned, we want the dimension of the boundary 
operator to remain fixed as we vary a and so we take m? = —?>/L 2 e . Thus A + = 3 
for our choice of mass, and we compute the solutions of ( P . 1 3| — p7T6|) numerically, 
reading off the r -3 fall-off of the scalar field to obtain (C 3 ) for a range of different 
temperatures. 
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3. Backreacting superconductors 



As shown in [|nj, the scalar field condenses out of its vacuum near the horizon 



of a sufficiently small black hole. Although the mass chosen for the scalar in [[15] 
is different to that used here (m' 2 = —3/L 2 > —3/L 2 ) we still expect a similar 
qualitative behaviour. Thus, in computing the dependence of (0 3 ) with T, we expect 
to see the characteristic curve depicting the condensation of (O3) from being trivially 
zero above some critical temperature T c , to some nonzero value below this critical 
temperature. Before proceeding with a full numerical analysis, it is useful to obtain 
some analytic understanding of the phase transition and the critical temperature. 

In order to find analytic bounds for the critical temperature, we look at the scalar 
field equation near T c . For temperatures just below T c , the scalar is only marginally 
away from its vacuum, hence the metric and gauge field will have the form ( j2.8| , 
|2~^ ), up to corrections of order 0(ip 2 ), thus the scalar field satisfies a linear equation 
( |2.14p with / and <ft taking their background values. 



A crude upper bound for T c can be found by considering the variable X = r 2 ijj, 
which satisfies to leading order: 



X''+[&-*)x'+(^ + -%--^)X = 0. (3.1 



fo rj V fo 2 LIU r/, 



At the horizon, 



which is positive for small k 2 Q 2 (taking X(r + ) > without loss of generality). Since 
A~l/rasr— >■ 00, the solution must have a maximum for some r, which requires 
that 

q 2 <\) 2 3 2f'\ 

+ ~m T > ( 3 - 3 ) 



f 2 L 2 J rhj 

at this point. An examination of when this condition is violated provides an upper 
bound for T c . This bound is only reliable for weakly gravitating systems, and figure 
H shows the upper bound for no backreaction, and k 2 = 0.05. For larger values of k 2 , 
the method fails because the behaviour of the function in ( p.3|) qualitatively changes. 

In order to obtain a lower bound, consider instead Y = rX, then manipulating 
the equation of motion for Y implies that if a solution exists, then the integral 



00 



« + i. + 3A_Ml = _rA^< (3.4) 

fo L 2 r 2 r \ >v 



r 3Y 2 



is negative. Note that negativity of this integral does not imply existence of a solution 
to the linearized equation near T c , it is simply a necessary condition. Since this 
integral is always negative at large T, and positive as T — > (for k 2 < 0.4), observing 
where it changes sign provides a lower bound on T c . This bound gives an extremely 
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reliable indicator of T c for a good range of k 2 , well within the numerical range we 
were able to explore. 

The analytic bounds for T c have been plotted in figure [I] together with the exact 
values of T c obtained by numerical computation. The upper bound has only been 
shown for k 2 < 0.05, as above this value it becomes less predictive and clutters the 
plot, and indeed beyond k 2 ~ 0.2 (corresponding to q ~ 2.25 in the notation of 0) 
it ceases to have quantitative value for any a. The lower bound on the other hand 
becomes successively more accurate as the values of h 2 are stepped up, and gives a 
very good quantitative guide to the behaviour of T c as we vary a and k 2 . 

It is easy to see that the effect of backreaction is to decrease T c and thus make 
condensation harder. We can see this from the point of view of the bounds, as 
backreaction at fixed Q drops the temperature of the black hole and also gives a 
smoother profile for fo, which makes it harder for the bound ( |3.4| ) to be satisfied. 
Essentially, the effect of backreaction is that the condensation of the scalar field not 
only screens the charge of the black hole, but also its mass, as the scalar and gauge 
fields now contribute to the ADM mass. This means that for a given charge and 
temperature, the radius of the black hole is increased, which makes it harder for the 
scalar to condense. 

One very interesting feature clearly exhibited in the bounds is the turning point 
in T c as a function of a for fixed k 2 . In the probe case, this occurs very near the Chern- 
Simons limit a = L 2 /4, and is barely perceptible in the numerical data, however this 
is shown clearly in the upper bound in particular. (We checked the numerical data 
by taking very small steps in the a parameter near L 2 /4, and were able to confirm 
this feature in the numerical data, however, these points are not shown for clarity.) 
Once backreaction is switched on, the minimum becomes much more pronounced, 
and indeed for large backreaction (k 2 = 0.2) the Chern-Simons limit is showing a 
considerable enhancement of T c over the typical values for lower a. 

In order to find the actual behaviour of the bulk superconductor, we integrated 
( |2 . 1 3| — pTT6|) numerically. As already stated, we took L — Q — q — 1, and varied 
r + to study how the system reacted to varying temperature. After solving for the 
matter and gravitational fields, the temperature is computed and the value of (O3) 
read off. 

Figure |^ shows (O3) 1 / 3 as a function of temperature for a variety of values of a 
and k 2 . Each line in the plot forms the characteristic curve of (O^) 1 ^ 3 condensing at 
some critical temperature. For simplicity we chose three values of a to display the 
features of the system: the Einstein limit (a = 0), the Chern-Simons limit (a = 0.25), 
and the mid-point value a = 0.125 to represent a 'generic' value of a. We also only 
present the detailed information for two values of k 2 : the case of no backreaction 
(k 2 = 0), and a backreaction of k 2 = 0.1. The plots are qualitatively similar for other 
values of k 2 and a, and can be deduced from the information in figure [I] on critical 
temperature. 
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Figure 1: A plot of the critical temperature as a function of a for a selection of k 2 . The 
analytic bounds are shown as lines and the numerical data as points. Respectively: k 2 = 
is shown in black, with solid lines and triangular data points; k 2 = 0.05 has blue dashed 
lines and diamonds; k 2 = 0.1 has a grey dotted line and squares; k 2 = 0.2 has a red dot- 
dash line with circular data points. The lower bound is shown for all n 2 values, but the 
upper bound is shown only for the lowest two values of k 2 , as they overlap significantly 
with the other data and confuse the plot. 

Figure @(a), which displays the raw data, shows how varying a and k 2 effects 
the height, shape and critical temperature of the condensate. This gives an alternate 
visualisation to figure p] of the effect of both backreaction and higher curvature terms 
on the critical temperature. Note how the backreacting case clearly exhibits the 
special nature of the Chern-Simons limit. Note also the disparity in the unnormalized 
condensate value. As in figure ^ we see that increasing k 2 reduces the critical 
temperature of the system markedly. We also see the effect of a, which initially 
lowers T c , then increases it again for a approaching the Chern-Simons limit. Because 
we are only showing the data for three values of a, the non-backreacting data do 
not display this effect, but it is clearly present in the backreacting data, and indeed 
increasing k 2 enhances this effect, 

Figure 0(b) shows the curves normalized by T c . The effect of k 2 is to increase 
the height of these graphs, in spite of the fact that the raw data tends to have lower 
values of (0 3 ). This is clearly because the most significant impact of increasing 
gravitational backreaction is that the critical temperature of the system is lowered. 
In this case the effect of backreaction is extremely marked, with the condensate 
varying more widely with a. 



4. Conductivity 



In |fl, Horowitz and Roberts observed an interesting phenomenon for the conductiv- 
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(a) (b) 

Figure 2: Two plots of the condensate as a function of temperature for a selection of 
values of a and k 2 . In each case, solid lines correspond to k 2 = and dotted lines to 
k 2 = 0.1. The black plot is a = 0, green is a = 0.125 and magenta is a = 0.25. The first 
plot shows unnormalized data, which indicates the variation of critical temperature as both 
a and k 2 vary. Plot (b) shows the conventional plot of condensate against temperature, 
both rendered dimensionless by normalizing to T c . 



ity of the boundary theory. They considered the model (|2.1| ) in the probe Einstein 
limit for a range of different bulk scalar masses, and on computing the conductivity 
found an apparent universal relation 

^ ^ 8 , (4.1) 

J- c 

with deviations of less than 8 %. In [|15| we found evidence that this "universality" 
was not in fact stable to the presence of stringy corrections, and we now wish to test 
the robustness of this result to backreaction. 

Conductivity is conventionally expressed as the current density response to an 
applied electric field: 

* = f (4-2) 

As the bulk field corresponds to a boundary four-current we must examine 
perturbations of to compute the conductivity. Since we are dealing with full grav- 
itational backreaction, we must also perturb the metric, and compute the variation 
of the EGB gravity equations. After some algebra we find: 

2- •■ L 2 fe 2u ( a(2u'f + f)\ 2K 2 r 2 A i 6' . . 

77 l r f e " A 'i}' + ^j AA - - ji 2 ^ A < + ~pr„ - l h " ~ ft-) = < 4 - 4 ) 

where h a b is the perturbation of the metric tensor, and Ai is the perturbation of the 
gauge field, which has only spatial components. Writing Ai(t, r, x l ) = 74(r)e ik ' x_;itJ 'ej, 
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and setting k = 0, we can integrate (1.3), and substitute in (4.4) to obtain: 

A = . (4.5) 



/ 



u; 2 __ 2 2 2 2/t 2 rV 2 
T 2 ^ 7 /e 2 " (r 2 - 2a f) 



We solve this under the physically imposed boundary condition of no outgoing 
radiation at the horizon: 



A(r) ~ f(r) 4 - T + . (4.6) 

Here, T + is the Hawking temperature defined at r = r+. In the asymptotic adS 
region (r — >• oo), the general solution takes the form 



where ao and a 2 are integration constants. Note there is an arbitrariness of scale in 
the logarithmic term, as pointed out in ||, however, this is related to an arbitrariness 
in the holographic renormalization process (see Appendix |A|) and we present here 
the expression used in our numerical computations to extract the behaviour of the 
gauge field. 

To calculate the conductivity, we must therefore compute the current dual to the 
gauge field fl4.7Q , and its linear response to an applied electric field. Since the details 
of the correspondence are a dependent (the asymptotic adS lengthscale changes 
with a), we go through the explicit calculation in an appendix. The result of this 
computation is 3 

Note that the imaginary term linear in u has an arbitrariness of scale from the 
counterterm subtraction, and in practise when we present the plots of u we will 
make use of this fact to choose an appropriate renormalization scale to give the 
greatest transparency of the features inherent in the conductivity. 

Figures || and f| present various aspects of our results for the conductivity at 
a range of values of a and k 2 . In figure |3| we show the real and imaginary parts 
of the conductivity as a function of oo/T c for no backreaction and a backreaction of 
k 2 = 0.05 for our three sample values of a: 0, 0.125, and 0.25. 

In each of these plots, the gap is clearly indicated by a rise in the real part of a, 
which coincides with the global minimum of Im(cx). As was noted in the derivation 
of ( |4.8j ), the imaginary part is only valid up to a linear term in u, the size of which 
is dependent on the renormalization scheme employed (and also on the charge Q). 
Certainly, within the parameter range tested, one can always find an appropriate 
term to create a finite global minimum in Im(cr) if it is not initially present. As 



3 Note that the equation for conductivity in jD| was missing some factors of L e : ( |4. 7| ) and (4 
correct equations (4.4) and (4.7) there. 
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in ||, we suggest that the value of this minimum be taken as u g , the value of the 
frequency gap. This has a clear advantage over trying to determine the midpoint 
of the step in Re(cr), in that it gives an unambiguous value of ui g . As the plots 
clearly show, once backreaction is included, the 'step' becomes much more gentle 
and extended, and while the dip in Im(cx) is smoothed, it is still clearly apparent, 
even for the most extreme, backreacting Chern-Simons, case. 




Figure 3: Conductivity: a range of plots showing the real (solid line) and imaginary 
(dashed line) parts of the conductivity as a function of frequency. Both variables are 
normalized by T c to dimensionless parameters. In each case the conductivity is shown for 
no backreaction in black, and for a backreacting parameter k 2 = 0.05 in magenta. Each 
plot represents a different representative value of a: (a) is a = 0, (b) is a = 0.125, and (c) 
is a = 0.25. The slight undulations in the plots at large w is a numerical artefact due to 
the sensitivity of the system near the horizon. 



The frequency gap is a distinct characteristic of a superconductor and in the 
BCS theory of superconductivity this frequency gap corresponds to the minimum 
energy required to break a Cooper pair. As mentioned above, in [|] it was claimed 
that for the holographic superconductor the relation oj g /T c ~ 8 had a certain uni- 
versality, proving stable for a range of scalar masses and dimensions. In this 
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relation was shown to be unstable to Gauss-Bonnet corrections in the probe limit 
(for m? = —3/L 2 ), and we unequivocally confirm this feature. Figure § gives a very 
clear indicator of how backreaction and higher curvature terms affect the gap. In- 
creasing either a or k 2 increases u g /T c . For the case of increasing a, the effect occurs 
mainly because of a shift in the gap, rather than a significant alteration of T c , which 
varies much more strongly with backreaction than a. On the other hand, varying k 2 
practically does not alter u g at all, whereas T c drops dramatically, leading to a sharp 
rise in u g /T c . 



i 




Figure 4: The gap frequency as a function of T c (the line uj g = 8T C is shown in black). Each 
data point on the graph represents a single pair (o,k 2 ). The different colours represent 
different degrees of backreaction; from right to left: Red is k 2 = 0, Blue is k? = 0.05, Green 
0.1, and Orange k 2 = 0.2. In each case a is incremented from to 0.25. As the gap alters 
rapidly near the Chern-Simons limit, the dotted lines are added by hand to guide the eye. 



5. Conclusion 

The aim of this paper was to understand how including both Gauss-Bonnet correc- 
tions and backreaction might affect the holographic superconductor. Our results are 
clear: increasing backreaction lowers the critical temperature of the superconductor 
hence increasing u g /T c . The effect of higher curvature terms is more subtle. Al- 
though these initially act in a similar fashion to backreaction in lowering the critical 
temperature, for significant GB coupling the critical temperature eventually begins 
to increase. The conductivity gap is also modified, with both ui g and T c altering to 
increase the ratio oo g /T c . We have therefore unambiguously refuted the claim that 
there is a universal gap u g ~ 8T C for these holographic superconductors, as even in 
the Einstein limit there is no such relation. 
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Our results show that there is a rich structure in higher dimensional holographic 
superconductors, with or without higher curvature corrections. For simplicity, we 
focussed on a single value of the scalar mass which fixed the dimension of the bound- 
ary operator. We believe this will clearly differentiate the effects of the backreaction 
and higher curvature coupling. From the results of [[L8j we do not expect any qual- 
itative differences to appear from a varying mass, although there will undoubtedly 
be quantitative differences, particularly since the expressions for conductivity used 
are inaccurate (as can easily be seen from dimensional grounds). It would be useful 
to explore more fully this parameter space, and to see if altering the bulk potential 
further enhances the features of the system, as well as exploring nonabelian gauge 
fields and more complex superconductors. 
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A. Holographic renormalization 



In this, we follow the method of Skenderis [f23| and explicitly compute the countert- 
erm and current dual to the electromagnetic bulk perturbation. First, we choose 
coordinates so that 4 

ds 2 = h^d^dC - L\% 

P 

Clearly, asymptotically, p = L 2 e /r 2 — > 0, and £ M = x M /L e with 



(A.l) 



ds 2 



dr 2 - dQ dp 2 



P 



P' 



(A.2) 



We then expand 



h 



A» = L-^ 2 [Af +P Af + ...] 



(A.3) 
(A.4) 



and solve the equations of motion order by order in p. (The factor of L e 1//2 ensures 
that the gauge fields Aft have the correct dimensionality.) 



*Note that the here is distinct from the metric perturbation notation used in section ^! 
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Focussing on the electromagnetic contribution, we can evaluate the action on- 



shell: 



s = [ -\f1^x 

Jm 4 

= / \A b V a F ab ^x - / \F ab n a A b Vhd< 

J M * JdM ^ 



X 



A'A vl ^ u ^)d A x 



(A.5) 



p=€ 



In these new coordinates, the gauge field equation of motion takes the form 

% ~ \^ Ai = ° {/) (A - 6) 

where d 2 ^ represents the wave operator with respect to the boundary metric 7^. 
The solution to this equation is 

A^L^AfHpL^ 

Thus the on-shell action, 



Af^-^pd^At 



(A.7) 



S 



A 



(0) 



^ 2) + ^)^ 0) + llne^' 



x 



(Ai 



is logarithmically divergent. In order to find the correct counterterm, we must invert 
the series solution ( |A.7|) to give A^ = y/L^Ai + O(e), and hence obtain: 

S rt = ^\ne [ Aidf^Ay^d^x 



p=e 



^lne / -FL Vhd 4 x 
4 7 p=e 2 



(A.9) 



We can now compute the boundary current, which is given by the 1-point func- 



tion 



5A, 



(A.10) 



Varying S + S ct explicitly gives 

5S + 5S ct = [ -F ab d a 5A b ^d 5 x + ^ lne f F^d^5A u Vhd 4 
Jm 2 J p=e 

= L e J ^d*x5A t (-24 + ^ df 0) JUj 



(A.11; 



Substituting for Aj from (|A.7| ) gives 



LI 1 / 



1/2 °S r , 



SAi 



A 



uiw f) 2 A {0) + -B 2 A {0) 



+ a (o) A i 



2 (°) * 



-2A< 2 > + 



ob 2 -k 2 



A 



(0) 



(A.12) 
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where the hatted quantities correspond to frequencies or wavenumbers with respect 
to the dimensionless coordinates £\ Notice that a shift in the renormalization scale 
e —7- Ae, results in a shift of the coefficient of the final term of In A/2. 

To get the conductivity, we can either compute the two-point function following 
p4}| , or use the standard formula (|4.2|) , noting that E = Af\ to get the dimensionless 



conductivity as 



2A& 



a 



100 



k=0 



+ Y 



(A.13) 



Finally, to obtain a dimensionally correct expression in terms of our original 
coordinates, we rewrite ( |A.7| ) 



At = L-^A^ + 



L 



3/2 r 



T 



Li 1 ' 2 



A 



(2) 



AP + S ( A® 



h^—B 2 4 (0) 
2 L (0) * 



(0) 



Lj(k 2 -u 2 
2r 2 



(A.14) 



ln(^)A (0) 



from which we obtain 



2a 2 . L e iu 
a = —— zwln(— ) + — 



(A.15) 



This is now the dimensionally correct expression for the conductivity, although there 
is an ambiguity in the imaginary part of iuj In A/2 as already noted. This differs from 



the expressions in [|T|, [L8| by factors of L e , mainly due to an incorrect extrapolation 
of the relation G R = — rfAA'lr^oo from the Einstein limit. 



References 

[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 
1113 (1999)] [arXiv:hep-th/9711200]. 

[2] C. P. Herzog, J. Phys. A 42, 343001 (2009) [arXiv:0904.1975 [hep-th]]. 

S. A. Hartnoll, Class. Quant. Grav. 26, 224002 (2009) [arXiv:0903.3246 [hep-th]]. 
G. T. Horowitz, "Introduction to Holographic Superconductors," arXiv: 1002. 1722 
[hep-th] . 

[3] J. D. Bekenstein, "Black hole hair: Twenty-five years after," arXiv:gr-qc/9605059. 
T. Hertog, Phys. Rev. D 74, 084008 (2006) [arXiv:gr-qc/0608075]. 

[4] P. Breitenlofmer and D. Z. Freedman, Annals Phys. 144, 249 (1982). 

[5] K. M. Lee, V. P. Nair and E. J. Weinberg, Phys. Rev. Lett. 68, 1100 (1992) 
[arXiv:hep-th/91 11045]. 

K. M. Lee, V. P. Nair and E. J. Weinberg, Phys. Rev. D 45, 2751 (1992) 
[arXiv:hep-th/9112008]. 



- 16 - 



[6] S. S. Gubser, Class. Quant. Grav. 22, 5121 (2005) [arXiv:hep-th/0505189]. 
Phys. Rev. D 78, 065034 (2008) [arXiv: 080 1.2977 [hep-th]]. 

[7] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, Phys. Rev. Lett. 101, 031601 
(2008) [arXiv:0803.3295 [hep-th]]. 

S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, JHEP 0812, 015 (2008) 
[arXiv:0810.1563 [hep-th]]. 

[8] G. T. Horowitz and M. M. Roberts, Phys. Rev. D 78, 126008 (2008) 
[arXiv:0810.1077 [hep-th]]. 

[9] S. S. Gubser and A. Nellore, Phys. Rev. D 80, 105007 (2009) [arXiv:0908.1972 
[hep-th]]. 

G. T. Horowitz and M. M. Roberts, JHEP 0911, 015 (2009) [arXiv:0908.3677 
[hep-th]]. 

R. A. Konoplya and A. Zhidenko, Phys. Lett. B 686, 199 (2010) [arXiv:0909.2138 
[hep-th]]. 

[10] S. S. Gubser and S. S. Pufu, JHEP 0811, 033 (2008) [arXiv:0805.2960 [hep-th]]. 
K. Peeters, J. Powell and M. Zamaklar, JHEP 0909, 101 (2009) [arXiv:0907.1508 
[hep-th]]. 

[11] K. Maeda, M. Natsuume and T. Okamura, Phys. Rev. D 79, 126004 (2009) 
[arXiv:0904.1914 [hep-th]]. 

[12] E. Nakano and W. Y. Wen, Phys. Rev. D 78, 046004 (2008) [arXiv:0804.3180 
[hep-th]]. 

T. Albash and C. V. Johnson, JHEP 0809, 121 (2008) [arXiv:0804.3466 [hep-th]]. 

[13] T. Takahashi and J. Soda, Phys. Rev. D 79, 104025 (2009) [arXiv:0902.2921 [gr-qc]]. 
S. Kanno and J. Soda, "Stability of Holographic Superconductors," arXiv: 1007.5002 
[hep-th]. 

[14] S. S. Gubser, C. P. Herzog, S. S. Pufu and T. Tesileanu, Phys. Rev. Lett. 103, 
141601 (2009) [[arXiv:0907.3510 [hep-th]]. 

J. P. Gauntlett, J. Sonner and T. Wiseman, Phys. Rev. Lett. 103, 151601 (2009) 
[arXiv:0907.3796 [hep-th]]. 

[15] R. Gregory, S. Kanno and J. Soda, JHEP 0910, 010 (2009) [arXiv:0907.3203 
[hep-th]]. 

[16] D. Lovelock, J. Math. Phys. 12, 498 (1971). 

[17] D. J. Gross and J. H. Sloan, Nucl. Phys. B 291, 41 (1987). 

R. R. Metsaev and A. A. Tseytlin, Nucl. Phys. B 293, 385 (1987). 

[18] Q. Pan, B. Wang, E. Papantonopoulos, J. Oliveira and A. B. Pavan, Phys. Rev. D 
81, 106007 (2010) [arXiv:0912.2475 [hep-th]]. 



-17- 



X. H. Ge, B. Wang, S. F. Wu and G. H. Yang, JHEP 1008, 108 (2010) 
[arXiv: 1002.4901 [hep-th]]. 

Q. Pan and B. Wang, "General holographic superconductor models with 
Gauss-Bonnet corrections," arXiv: 1005.4743 [hep-th]. 

R. G. Cai, Z. Y. Nie and H. Q. Zhang, "Holographic p-wave superconductors from 
Gauss-Bonnet gravity," arXiv:1007.3321 [hep-th]. 

[19] Y. Brihaye and B. Hartmann, Phys. Rev. D 81, 126008 (2010) [arXiv: 1003.5130 
[hep-th]]. 

[20] D. G. Boulware and S. Descr, Phys. Rev. Lett. 55, 2656 (1985). 
R. G. Cai, Phys. Rev. D 65, 084014 (2002) [arXiv:hep-th/0109133]. 

[21] L. F. Abbott and S. Deser, Nucl. Phys. B 195, 76 (1982). 

S. Deser and B. Tekin, Phys. Rev. D 67, 084009 (2003) [arXiv:hep-th/0212292]. 

[22] M. Banados, C. Teitelboim and J. Zanelli, Phys. Rev. D 49, 975 (1994) 
[arXiv:gr-qc/9307033]. 

[23] K. Skenderis, Class. Quant. Grav. 19, 5849 (2002) [arXiv:hep-th/0209067]. 

[24] D. T. Son and A. O. Starinets, JHEP 0209, 042 (2002) [arXiv:hep-th/0205051]. 



- 18 - 



